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ABSTRACT 


Under laboratory conditions particle losses from a magnetized plasma will produce a neutral 
particle flux reentering from the surrounding vessel walls into the plasma region. Even at very 
high temperatures, where the plasma is almost fully ionized, this gives rise to interactions between 
charged and neutral particles which cannot always be neglected. Three types of interactions are 
discussed; ionization, charge exchange and elastic collisions with neutrals. 

The conservation laws of mass, charge, momentum and energy are reconsidered for a two- 
fluid system with ionization, charge exchange and elastic scattering of neutral particles taken 
into account. It is found not to be necessary to derive the energy equation in terms of the joint 
action of electromagnetic and mechanical changes of state. 

It has earlier been concluded from the motion of a single particle that a very strong confinement 
can be obtained in a rotating plasma situated in the magnetic field of a current loop (Bonnevier 
and Lehnert). These results are verified in a macroscopic theory including the energy balance of 
stationary motion. It is found that the main part of the heating of the plasma and the scattering 
of ions are not provided by coulomb collisions, but by collisions with reentering neutral particles. 
Two special cases are treated in detail, one with vanishing heat conductivity and one with infinite 
heat conductivity in the direction along the magnetic field. For the latter, which is more likely 
to occur in a high temperature plasma, the energy loss by escaping particles is deduced. The 
theory predicts that this loss should be several orders of magnitude less in the current loop con- 
figuration than in earlier described homopolar machines. 

Finally, some comments are given on the non-existence of interchange instabilities in systems 


with volume currents. 


I. Introduction 


The ultimate goal of controlled fusion research is to produce and confine a plasma 
at extremely high temperatures. For this purpose the flux of high energy particles 
escaping from the plasma to the surrounding vessel walls should be kept at the lowest 
possible level. Otherwise the thermal coupling between the plasma and the walls 
becomes very large and, drastically speaking, the plasma does not reach the intended 
high temperature unless the vessel walls would be “boiled away’”’. 

As an example, there is only hope to balance the energy loss of escaping particles 
at some 10° °K by the reaction energy in a deuterium plasma with density n = 10°! m 
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if the particles are confined during times of the order of a second. Ions then have to 
travel about 10° m before they leave the confinement region. Without any confine- 
ment the energy flux to the vessel walls would become about 10! watts/m?. Suppose 
that the confinement is produced by a magnetic field with a mirror ratio of about 2. 
The velocity distribution of the particles then becomes anisotropic. Taking only 
coulomb collisions into account the confinement time becomes about 10~* sec. and 
the energy loss by escaping particles about 10° watts/m* (Post 1958). At 10° °K the 
corresponding values become 10~° sec. and 10° watts/m*. Only if temperatures around 
3 x 10° °K are reached, a mirror machine may become self-sustained (Post 1958, 
Garren et al. 1958). 

Further complications in the energy balance of a thermonuclear fusion device are 
introduced by such effects as synchrotron radiation (Trubnikov and Kudryatsev 
1958), interactions with neutral gas entering the plasma and by instabilities. 

In the present paper the energy balance is studied for a system with small particle 
losses. An application is‘given by the plasma confined by the magnetic field of a 
current loop as suggested by Lehnert (1958a) and later by Greyber (1958) and 
Colgate and Furth (1958), heating of the plasma being provided by a transverse 
electric field which sets the plasma into rotation around the axis of symmetry (Leh- 
nert 1959). In an earlier study of the motion of an individual particle it has been 
shown that particles under certain conditions become trapped completely inside a 
“bag” of magnetic field lines which is ‘‘sealed’”’ by the centrifugal force (Bonnevier 
and Lehnert 1959). The situation resembles that of a planetary atmosphere, where 
the particles are trapped by a gravitation field, provided that their thermal velocity 
falls below the escape velocity. To determine the situation completely the relation 
between the thermal energy of the plasma and the centrifugal force has to be deduced 
from the energy balance. A macroscopic theory is given for this purpose in the present 
paper which also presents a determination of the energy losses by particles escaping 
from the plasma. 

Some astrophysical applications to the problem of a rotating plasma are discussed 
more fully elsewhere (Lehnert 1960). 


II. Conservation laws for the two-fluid model 


Losses of charged particles from a magnetized plasma to the walls of a surrounding 
vessel will produce a flux of neutral particles reentering the plasma region (see also 
Luce 1958 and Miller 1949). Even if the plasma temperature is very high this implies 
that the system behaves like a partially ionized gas. Conservation of momentum for 
such a gas has earlier been considered by Schliiter (1950, 1951), Lucas and Schliiter — 
(1954), Cowling (1956, 1957), Piddington (1956), Gershman and Ginzburg (1956) and 


Lehnert (1958 6). In the present paper conservation of energy will be included in the 
treatment. 


1. Collisions between charged and neutral particles 


In this paper the conservation theorems will be extended to include ionization 
processes and charge exchange collisions. Throughout the paper will be assumed that 
the charged particles which escape to the walls will recombine there and reenter into 
the gas volume in the form of neutral particles with a thermal velocity distribution 
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nearly corresponding to the wall temperature. The possibility of using a wall material 
with high adsorption ability will not be discussed here. 

Further, the treatment is simplified such as not to include the spontaneous emis- 
sion of photons by excitation, bremsstrahlung and synchrotron radiation as well as 
two-stage ionization, volume recombination and multiple charge exchange collisions. 
The wall temperature is assumed to be considerably lower than that of the plasma. 
Finally, the plasma is assumed to be electrically quasi-neutral with n ions and elec- 
trons per unit volume. 

Under these restrictions the reentering neutral particles may react with the plasma 
in three different ways. Firstly, ionization may take place with a number of ionizing 
collisions per unit volume and time 


NM Ny (0; Uj» +N Mn (Oe Ue =N* E, (1) 


where 1, is the local neutral gas density inside the plasma, o; and o, the cross sections 
for ionization by means of ions and electrons, u; and u, the total velocities of ions 
and electrons, ¢ > indicates effective mean values and é expresses the number of charged 
particles created per unit volume and time per charged particle being present. The 
second process to be taken into account is that of charge exchange between ions and 
neutral atoms. With the cross section given by o,, and the corresponding collision 
frequency by y,, the number of charge exchange collisions per unit volume and 
time becomes 


NMNn (Cer UW) =N* Ve, = N° f &, (2) 
where eq. (2) is to be taken as a definition of f. Thirdly, elastic collisions with cross 
sections o;, and o,, may take place between charged and neutral particles with rates 
given by 

NM Nn Fin Uj» =20 Vin =22 Gin & (3) 


and n Nn <Oon Ue» ait Ven <a NGJen & (4) 


where »,, and y,,, define the effective collision frequencies between ions and neutrals 
and electrons and neutrals, respectively. In average ions lose half of their momentum 
relative to the neutrals. The mean free path for ionization of a neutral particle which 
enters the plasma region with the velocity u,, becomes 


DL, = My Uyn|n &, (5) 

and the corresponding paths for exchange collisions and elastic collisions become 
Lys =. i / f (6) 

and Lin = Li | Gini Den = Li / Gen: (7) 
Finally, we also introduce the rate at which neutral particles are scattered elasti- 

ca lly by ions, thereby obtaining a mean energy comparable to that of the ions and being 
lost to the vessel walls before undergoing any ionizing or charge exchange collision. 
The number of such processes per unit volume and time will be indicated by n-h & 


where A nearly equals g;,/2. A process of this kind may take place when L,, is 
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comparable to L, and L,;. Then, the neutral particle may be “kicked” out of the 
plasma before it has been involved in any ionizing or charge exchange collision. 
This provides an additional energy loss mechanism with the same net result as 
that of the charge exchange. In connexion with the “DCX” experiments Barnett 
et al. (1958) have discussed the conditions for “burning out”’ the neutral gas in- 
side a plasma, where neutral particles are subject to ionization and charge exchange. 
In cases where the mechanism just mentioned comes into operation the conditions 
for “burnout’’ may have to be reconsidered. 

It should be pointed out that, before being ionized, part of the reentering 
neutral gas flux may become subject to elastic collisions with electrons. As a 
result the neutrals will diffuse through the electron gas, thereby being slowly 
heated and accelerated until ionization or collisions with ions take place. A 
detailed analysis of this situation will not be undertaken here. 

Thermal diffusion effects are not included in the present treatment. 


2. Conservation of mass and charge 


With the starting points given in the preceding section the conservation laws will 
now be formulated for a quasi-neutral plasma. The mass density of the ion gas is 
given by 0; =nm,, the temperature by 7';, the pressure which is assumed to be iso- 
tropic by p; =nkT', and the mass velocity by v,. Corresponding quantities for the 
electron gas are indicated by the subscript (,). Also introduce the total plasma 
density 

o=n(m,+m,) =nM = 0; + Oe, (8) 
the plasma pressure 
P~ Pi Ze Pe> (9) 


the mass velocity v of the plasma and the electric current density i. The two latter 
quantities obey the relations 


Vv; =V+m,i/e0; V.=V — m,i/eo. (10) 


Since charged particles are created only by ionization conservation of mass and 
charge is expressed by 


div (wv,) =En SE = div (nY.). (11) 


Introduce the derivative d/dt =0/ét +v-Y and combination of eqs. (10) and (11) 
gives 
: 1 do 
div v= &-—--— 
iv nee (12) 


and div i=0. . (13) 


Consequently, the starting points also require that the displacement current has to 
be neglected and 


curl B = py i, (14) 
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where 4) is the permeability in vacuum. This is a good approximation when the 
characteristic velocities of the phenomena in question are much less than the velocity 
of light. The diamagnetism of the plasma can be taken into account without introduc- 
ing a modified permeability; the current density i in eq. (14) expresses all currents 
within the plasma explicitly (Spitzer 1956). 


3. Conservation of momentum 
Throughout this paper will be assumed that the major part of the neutral gas 
which interacts with the charged particles has a negligible mass velocity and a low 
temperature with respect to the plasma. These conditions will be discussed further 


in section I1.5. The gain of momentum by motion in the electric and magnetic fields 
E and B then becomes 


i 


ra =en (E+ v; xB) —V Di —~ Oe Vei (Vv; =¥z) — Oi (g St Perce Vin) Vi (15) 
for ions and 
de 
Rea, Ve= —en(E+v.XB)—V pe+ Oe Ver (Vi — Ve) alot.) Ve (16) 
for electrons, where »,; is the frequency of coulomb collisions and 


d,/dt =d/dt+v,-V; d./dt=d/ét+v.-V. 


By making the substitutions (8)-(10) expressions (15) and (16) can be written 
in terms of vy and i as shown in an earlier paper (Lehnert 19580). The sum of eqs. (15) 
and (16) gives the conservation law of the mean motion of the plasma 


9 i iY) i/o)=ixB—V p—aov+Bi, (17) 
where a= (ME+M, Veg + M, Vin +M, Yen) |™M (18) 
and B =™M™; Me (Ven Tass Vex) (em (19) 


are “‘frictional coefficients” related to the loss of momentum by collision processes. 
If instead eq. (15) is multiplied by m,/eo and eq. (16) by — m;/eo the sum of the 
resulting equations becomes 


e d F ° mMm—™Me . . 
fae lee) + Gv) y— ; ‘(i-V) ‘/a) = 
F (20) 
- i] } 
=E+vxB— ix B+ — (m, V Be— me V 1) + BV—Ni, 
eo eo 
2 
MMe m M; Me m 
where qa apt eat e+ 2 (Yin + Vex) + (“) ra (21) 
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Kq. (20) is a generalization of Ohm’s law and expresses the conservation of momentum 
of the “slip” motion between the ions and electrons; 7 is a generalized electric re- 
sistivity. For further discussions of the significance of the different terms in eqs. 
(17) and (20) reference is made to earlier investigations (Lehnert 19585). 


4. Conservation of energy 


If scalar multiplication of eq. (17) by v and of eq. (20) by i is performed the sum 
of the resulting equations can be written 


mM Pe Me Pi 


E-i-— div (py Pe ) = ta+ w+ (22) 
Q 


The first term in the left hand member of eq. (22) represents the electric power input 
per unit volume. The second term in the same member is equal to — div (p;V; + .V-) 
and represents the total work done by the pressure on a fluid element. Thus, the left 
hand member represents the total energy input. The right hand member tells how 
this energy is consumed. Here 


d; de dy ™,mM... p 
Wa = OMIT > YEP CGY ely PY ores a ras e (i- V) (i/o) cis 


t 
(23) 


Mm; Me 


a be 1 Me GW) (i 
+ ot i fo S/o) + vy GW) G/a) 


is the total acceleration work, 
W. = — p; div V;— p, div Vv. = — p div v + (m;p, — m,‘p;) div (i/e @) (24) 


the compression work and 


where Wry = o(& a Ver "5 Vin) Vi aE Qe Vei(Vi 52) °v; (26) 
Wye ra 02(& rz Ven) Ve —~ @e Voi (Vj tn me be (27) 


expresses the frictional work done on ions and electrons, respectively. 

It should be observed that eq. (22) can be rewritten in such a way that the pressure 
p is everywhere replaced by p* = p + B?/2 jy, where the second term is the “mag- 
netic pressure’’. The right hand side of eq. (22) will then contain a “magnetic compres- 
sion work” —(B?/2 9) div v, an acceleration work —v-Y(B?/2 4) done by the 
magnetic pressure and an additional acceleration work v-[(B-V) B/ Mo] which is due 
to the fact that the magnetic field lines also may act like elastic strings. The force _ 
(B-V) B/uo usually does not have the character of a pressure gradient. However, this 
way of describing the energy balance is unnecessarily complicated and does not give 
any result beyond that of eq. (22), which is written in terms of the real pressures p; 
and p, and the total Lorentz force i x B. 
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The obtained relation does not give any new information beyond the laws of con- 
servation of momentum from which it is deduced. Information is still required about 
the coupling between thermal motion and mass motion and this can be obtained from 
the principle of conservation of energy (see also Lehnert 1956, 1958c, Pai 1957, 
Bafios 1958 and Chu 1959). 

First consider the energy balance separately for the ion gas inside a volume element 
which is fixed in space. The increase in total energy of the gas is 


7) 
ay le Pit 2: Vi)- 
The energy input becomes 

enk-vy; ae, div (p;V;). 


The plasma temperature will be assumed high enough to make the ionization work 
negligible as well as the neutral gas temperature. 
Convection will produce a power loss 


div (3p, Vv; +30, Vi- Vi) + $pi+40; Vi) (Yer thé), 


where the first term is due to motion of the ion gas and the second to energetic 
neutral particles escaping from the plasma region, / being defined at the end of sec- 
tion II.1. Heat conduction described by the flow vector Q; produces a power loss 
div Q,. Heat radiation will not be considered here. Finally, the ions will lose the power 
w,,; due to mechanical work performed against the frictional forces. The fraction 
0; of w,; is thereby used to heat the electron gas and the parts of the neutral gas 
which do not become subject to processes defined by »,, and A. Conservation of 
energy requires that 


pp lePit Qi Vi) =en Evi — div (p; V;) — div (39, %4 +40: Vi * Vi) — 
: 28 
~ 3p; + 40; V2) (ver th E)— div Q)— 4,09. CS) 


It should be stressed that e(E +u, x B)-u,; =eE-u, is the total energy input from 

the electromagnetic field into the motion of an ion which has the total velocity u,. 

Tf the magnetic field B = curl A changes in time this will produce an energy input 

by means of the induced electric field part —é A /ét. Thus, enE-v; is the only link by 

which energy can be fed into the ion gas by the electromagnetic field, irrespective 
of the kind of electromagnetic processes being in operation (Lehnert 1956). For the 
electron gas the corresponding relation becomes 


F ve+}o0v8)- —enE-Ve— div (pe Ve) — div (3 PeVe+}0cVe*Ve)— div Q.-dewye. (29) 


Eqs. (28) and (29) express the conservation of the total energy as seen from the point 
of view of the charged particles. The equations express that the net change in the 
thermal and kinetic energy of the plasma is caused by an outflow of particles, by 
heat conduction, and finally by energy which is transformed into other forms than 
thermal and mechanical energy. In this case the latter form of energy is electro- 
magnetic and the first term of the right hand sides in eqs. (28) and (29) gives the rate 
at which this energy is supplied from the electromagnetic field. “Radiation” of mag- 
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neto-hydrodynamic waves is also included in eqs. (28) and (29) since these waves are 
characterized by the work done on the charged particles by the induced electric field. 
Consequently, it is not necessary, at least in the case of an ionized gas, to deduce the 
energy equation in terms of the joint action of electromagnetic and mechanical 
changes of state as suggested by Chu (1959). 

The coupling with the electromagnetic field is now eliminated by means of the con- 
ditions for conservation of momentum. After scalar multiplication by v,; and v, eqs. 
(15) and (16) are subtracted from eqs. (28) and (29) and div v, and div v, are substitu- 
ted from eq. (11) into the obtained relations. After some deductions the result 
becomes 


d; (5 4 Pi Ui Qi _ 


i 4 +3 eat eh i+ Oe Ver (Vi — Ve) * Vi + 
di \2 a) O ode (ETB ver +305) Bi + Qe Pei (Vi Ve) Vi an 
+30; (E+ Veg—hE+2 Vin) Vi — div Q; — 6; wx; 
de (3 Pe\ Pe te Ge 5 
> ne ; are Fr le — Oe Vey (Vi — Ve) * Ve + 
gt (5 ) oe dt g & Pe Ge Pei (Vi — Ve) -¥ 
(31) 


+3 Oe (E+ 2 Yen) v-— div Q. — be Wye. 


If energy losses, heat conduction and the neutral gas inflow disappear eqs. (30) and 
(31) reduce to the simple adiabatic relation where the pressure is proportional to the 
density raised to the ratio between the specific heats. By means of the substitutions 
(8)-(10) the sum of eqs. (30) and (31) is transformed into the energy theorem 


d (3 3 WY TOLL Vd} 
(1—d)uy~ div @-2Ep— 8+ WED = [os (52) —Si-w (PMB) 


do m gi Lilla) Dir. 2 
- |P.58— ea "G-Ve)-Ep-L+/+ME ove -bEoer?|, 


where w, is the total heat power created by friction as given by eq. (25) and Q =Q, + 

Q), is the total heat flow by conduction. The left hand member is the net flow of heat, 
where the heat loss due to escaping charged particles is given by the third term and 
the heat loss due to energetic neutrals escaping from the plasma is indicated by the 
fourth term. The first bracket of the right hand side gives the rate of growth of the 
internal energy and the second bracket is the total mechanical work. The first three 
terms inside this bracket represent the compression work; £ p is the work which has 
to be done to “force’’ n& new particles into the plasma volume at the pressure p. 

The last two terms within the bracket give the mechanical work done by frictional’ 
forces to accelerate né new particles to the full mass velocity and to balance the 
momentum loss of the mass motion which escapes with neutrals out of the plasma 
volume. The heat lost from the plasma to the neutral gas, which does not escape 
directly to the walls with energetic neutrals, is given by 6 wy. 


Finally, introduce the total flow of energy q caused by convection and heat conduc- 
tion by means of charged particles; 


4 =3 (Di Vi + De Ve) +4 (0; 07 Vi + Oc V2 Ve) + Q. (33) 


The energy theorem (32) may then be rewritten by taking the sum of eqs. (28) 
and (29) and combine it with eqs. (22) and (11). The result becomes 
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Ai sPpt+rsov Reis = Wat Wetn(1 —d) w,— div q—(§7:+50, Vi): (f+h)é. (34) 


This relation describes that the growth of total energy given by the left hand member 
is balanced by the sum of the acceleration and compression works and the heating 
by friction reduced by the energy losses. The latter are produced by convection and 
heat conduction by charged particles together with energy losses caused by escaping 
energetic neutrals produced by charge exchange and by elastic scattering against 
ions. The present results are easily generalized to include a mass velocity of the 
neutral gas. 


5. Applications to gas discharge experiments 


When the theory is applied to gas discharge experiments on laboratory scale a 
number of approximations can be made. The following questions arise: 


a. “Creeping diffusion’’ and inertia forces acting on the electric current 


In the equation of motion (17) the second term of the left hand side represents the 
difference in acceleration of the ion and electron gases. When this effect is neglected 
the plasma is said to be in a state of “‘creeping diffusion” (Schliiter 1951, 1956). 
In Ohm’s law (20) the left hand side represents the inertia of the charges which 
produce the current. A rough estimation of orders of magnitude is now made by 
replacing space and time derivatives by inverted characteristic lengths and times. 
Throughout this paper the characteristic dimensions over which all quantities vary 
in space are assumed to be comparable with the dimensions of the confinement region. 
Thus, surface currents with infinite current density and pressure jumps across 
discontinuity surfaces are excluded. Then, it is easily seen that the second term in 
the left hand member of eq. (17) as well as the left hand member of eq. (20) can be 
neglected, provided that the time scale of the macroscopic phenomenon to be con- 
sidered is assumed to be much longer than the gyro time of an electron in the magnetic 
field B (cf. Lehnert 19585). 


b. Dissipation 


A strong confining magnetic field will be assumed to exist in the systems to be 
treated and the dissipation of energy as well as the particle losses and the influx of 
neutral particles are assumed to be small. A rough estimation also shows that the 
terms with coefficients § and 7 in eqs. (17) and (20) will become small when the fre- 
quencies £, ¥,;, v;, and ¥,, all are much less than the gyro frequency of an electron. 


ce. Collisions with neutral particles 


In cases where / defined in section II.1 is much less than unity a neutral particle 
entering the plasma region will have a good chance to be subject to ionization be- 
fore it gains momentum in an elastic collision with an ion and is scattered out of 
the plasma. The momentum gained by elastic collisions with electrons is about 
m,|m times less and can be neglected even if the velocity of the electrons is about 
(m,/m,)* times higher than that of the ions, provided that 7; and 7, do not differ 
too much. If at the same time at least one of the free paths LJ, and L,, defined in 
eqs. (5) and (6) is smaller than the thickness of the plasma region all the entering 
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neutral gas flux will be “absorbed” in the plasma and subject to ionization or charge 
exchange. This implies that the energy loss defined by h and 6 and the momentum 
loss given by v;, and »,,, disappear in all eqs. (17), (20), (25), (32)and (34); all the losses 
are already represented by terms including é and »,,. This simplified situation will 
be studied henceforth. 

The assumptions seem to be valid for deuterium in a temperature range from 
4x 10° to about 3x 10° °K. According to Fite (1958) the values of o, vary in this range 
approximately from 7 x 10-71 to 3x 10-?! m?. Since oj, 5.9 x 10-*° m2 the quantity 
h 1 Gin U;>|<Oo Ue> Will vary between the corresponding limits 0.035 and 0.090 when 
T,-~T,. Ionization by ions, which has not been taken into account, somewhat re- 
duces these values of h. It should be observed that in this temperature range the 
ionization rate is mainly provided by electrons and that it exceeds the charge ex- 
change rate. According to Little (1956) o,,<1.5 x 10-*° m? for 7’; in the range men- 
tioned above giving f < <o.,U;> /(d.Ue> < 9.10. 

Finally, a plasma density n> 1018 m-3 gives paths for ionization L,<0.05 m 
when 7',=2 x 10° °K and the neutral particles enter the plasma region with a 
velocity wu, corresponding to a wall temperature 7',, = 300 °K. In such a case it seems 
likely that the particles will be absorbed during a single passage through the plasma. 


III. Confinement and energy balance of a rotating plasma 


Consider a mass of ionized gas which is rotating around an axis of symmetry. The 
rotation is produced by an electric field E the part E, of which is perpendicular to the 
magnetic field B, as in the ‘““Homopolar’’ and ‘‘Ixion” devices discussed by Baker 
and Anderson (1956) and Boyer et al. (1958). The part of B generated by sources out- 


current lead 


Fig. 1. The shaded area in the figure represents the cross section of a thin toroidal shell generated 

by the field lines of a strong poloidal magnetic field B which is symmetric around the z-axis of 

the figure. An electric field E, essentially at right angles to B, produces a rotation around the 
z-axis for charged particles inside the shell (cf. Bonnevier and Lehnert 1959). 
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side the plasma is supposed to be purely poloidal. Introduce a cylindrical coordinate 
system (7, p, z) with the coordinate z along the axis of symmetry as shown in Fig. 1. 
Stationary, axially symmetric situations will be studied, where all macroscopic 
quantities are functions which do not depend explicitly on and t. The electric field 


becomes 
= 0¢ , @¢ 
E= / or’ 0, az; ‘ (35) 


The shaded area in Fig. 1 represents the cross section of a thin toroidal shell generated 
by the field lines of a strong poloidal magnetic field B. Confinement of charged parti- 
cles inside the shell will now be discussed both in terms of single particle motion and 
from the macroscopic point of view. 


1. Motion of an individual particle 


In a study of the forbidden regions for a single charged particle moving in the 
magnetic field of a current loop (cf. Fig. 1) Bonnevier and Lehnert (1959) have found 
that complete particle confinement can be achieved under certain conditions. This 
occurs when a transverse electric field is applied which produces a drift motion and 
a rotation around the axis of symmetry. As a consequence of the shape of the mag- 
netic field the centrifugal work becomes very large for a particle which is “‘lifted”’ all 
the way from the “‘starting point”’ (79, 0) in Fig. 1 to a point (r,, 0) in the region of the 
current leads which suspend the loop. 

The magnetic field of Fig. 1 is now assumed to be a vacuum field and an electric 
field is applied across the narrow space between two electrodes which follow two 
toroidal surfaces generated by the field lines. E and B = curl A then become orthogonal 
all over the space between the electrodes. The shaded strip of the figure contains the 


magnetic flux 
d® =2arBdl=2zxd(rA,) (36) 


and the potential difference between its inner and outer surfaces is 
d®= —-Edl. (37) 
The angular velocity of rotation inside the strip becomes 
Q= — E/Br=22dd/d® =dd/d(rA,). (38) 


Since d¢ and d® are constant along the strip this shows that Q is constant along a 
magnetic field line. With /u,, uw, +Qr, u,/ as the total velocity of a charged particle 
of mass m, and charge q, the relation 


uw + uz = ud — (12 — 7?) Q? — 2 (q,/m,) (6 — bo) + 2 Q (G/m,) (7 Ag— 1% Ago) — nn 
— [19 Upo— (qo/My) (Ag — 79 Ago) + (79 — 7") QP/r* 


has been obtained by Bonnevier and Lehnert (1959). For very strong magnetic fields 
the particle will be guided along a very narrow toroidal region like that indicated 
by the shaded area in Fig. 1. Relation (38) then shows that the third and fourth 
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terms in the right hand member of eq. (39) cancel. Then, the particle confinement 
becomes complete when the centrifugal work (m, /2) (76 — 12) Q2 exceeds or becomes 
equal to the particle energy (m,/ 2)us at the starting point (74, 0); the point r be- 
comes a turning point. 

In the case of a plasma there are interactions between particles which have not 
been taken into account in the treatment of the individual particle motion. Further, 
the degree of confinement produced by the centrifugal field cannot be determined 
unless a connection between the rotation and the thermal motion is established. This 
requires a study of the energy balance as shown by the next paragraph. 


2. Macroscopic theory 


When the plasma confinement is studied macroscopically for configurations of the 
type given in Fig. 1 the assumption of isotropic pressures p; and p, gives a situation 
which is less favourable than the real one, where particles are scattered into the loss 
cones by collisions. Thus, the energy losses calculated in this section are to be con- 
sidered as an upper limit. 

The problem to be solved is that of a plasma “atmosphere” contained in a “bag” 
of magnetic field lines as indicated in Fig. 1 and being rotated around the axis of 
symmetry. Calculation will be made of the rate at which this atmosphere escapes 
out of the “holes” of the magnetic “‘bag’’ which are here, for the sake of simplicity, 
represented by the somewhat less favourable situation of a ringshaped wall situated 
at r=r,. 


a. Basic considerations 


With the assumptions of section I1.5 of a system with small losses scalar multiplica- 
tion of Ohm’s law (20) by B gives a longitudinal electric field E,,, where 


E, Be —B-(mV Pe — mM. V pi) eo. (40) 


Consequently, the difference in electric potential energy along a field line and along 
the whole configuration becomes of the order of the thermal energy of the particles. 
It is now supposed that a strong electric field part E, is applied in such a way that 
the transverse potential difference across the configuration is much greater than that 
corresponding to the thermal energy of the particles. This is possible if a sufficiently 
strong magnetic field is applied to make the radius of gyration of a charged particle 
much smaller than the dimensions of the total confinement region. Thus, the situation 
to be studied is characterized by H,>E,, where the electric field is very nearly 
perpendicular to the magnetic field almost all over the plasma region with the 
exception of the Debye sheets. 
Combination of eqs. (17) and (20) gives 


m,—Me. Av 
E+vxBe + + (ms VD — eV Be)/e @ + 


+ [n — (m, — me) B/e e]i—[B — (m;— m,) «/e] v. 


A rough estimation of orders of magnitude immediately shows that the right hand 
member becomes much smaller than the second term of the left hand member when 
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the gyro times of the charged particles are much shorter than the characteristic time 
changes of the macroscopic motion and the collision times. This is a fair approxima- 
tion at small losses and strong magnetic fields. 

The assumption of small energy losses also implies that the energy input E-i 
should be small. Since H,=0 and £ is relatively large the current components i, 
and 7, must be small quantities, as well as the toroidal magnetic field b, generated 
by these currents. On the other hand, no restriction is made upon the induced toroidal 
current i, which gives an electromagnetic force that balances the centrifugal force 
and the pressure gradient in the equation of motion (17). The consequence of this is 
that the velocity of rotation becomes 


Up =Qr wz (E x B),/B. (42) 


Since E is almost perpendicular to B eqs. (36)-(38) show that Q becomes nearly 
constant along a field line. This is consistent with the conclusion earlier drawn 
by Ferraro (1937) about the isorotation of a magnetized star. 
Introduce the part v,,, of the transverse velocity v, which is situated in the rz- 
planes: 
Virz=V, —Vo BX (E+Vv x B)/B?. (43) 


From eqs. (41) and (43) is then immediately seen that v,,, can be made very small 
when the magnetic field becomes strong enough. This only implies that the transverse 
motion of matter across the magnetic field in the rz-planes becomes very small and 
that nearly all the flux of charged particles out of the confinement region may take 
place along the magnetic field lines with the velocity v,,, where v)>v,,,. Conse- 
quently, the equation (12) of continuity has the approximate form 


@:€ = div(ov) ~div(ov,)) ~ + div (Box, /B) = + B-V(0%,/B), (44) 


since div B = 0 and the plus sign is chosen when y,, is directed in the positive direction 


of B. 
The coordinate s along a magnetic field line is introduced for the upper half plane 


of Fig. 1 where 
a Bb, a earn B 


=—. eames: 45 
Tae ree at TD (49) 
Eq. (44) is integrated to 
$ ue 
Beier ae POs (46) 
V0); Bien | = dr=—TI(s); (20), 
i ab oJ B, Q 


where s =O corresponds to r=7y, which gives v,,=0 by reasons of symmetry. 
Integration is made along a field line. From this expression and eq. (42) is seen 
that v,, becomes much less than v, when the particle losses are small enough to make 
the necessary frequency & of ionization per charged particle much smaller than the 
angular frequency Q. This is the situation to be studied henceforth. It implies that 
a charged particle travels many times around the axis of symmetry before it is lost 
out of the confinement volume. 
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If the equation of motion (17) is multiplied scalarly by B the result for small 
losses becomes 
dp 


B-Yp=B-- ~ —eB-[(v-V) v]~oB, Q's, (47) 


i.e., the pressure gradient in the “‘atmosphere’”’ of the magnetic field “bag” is balanced 
by the centrifugal field. Especially when p; = p. = p/2 the longitudinal electric field 
of this atmosphere can be deduced from eqs. (17) and (41) which give 


E-Bw — (m, — m,) B, Q?r/2e. 48) 


Due to their smaller mass electrons have a tendency to “evaporate”? more easily 
from the atmosphere than the ions. The electric field E,, of eq. (48) is directed in a 
way such as to retard electrons and accelerate ions, thereby counteracting a charge 
separation and producing equal escape rates of both kinds of particles. A similar 
result has been obtained by Spitzer (1956) for a plane atmosphere balanced by the 
action of gravity. 

The part of the electric current which is perpendicular to B is obtained from eq. 
(17) and is given by 


Bi, ~o9 Bx[(v-V) V]+BxyV pe /° — 9 Q?rB,2P— pcb, 0 / (49) 
in the first approximation. Since 7, and 7, are small an estimation of orders of magni- 
tude by means of eqs. (42) and (49) shows that 7/env is of the order of the ratio 
between the radius of gyration of an ion and the dimensions of the confinement 
region, provided that thermal and rotation velocities are of the same order of magni- 
tude. An examination of the frictional work can now be made by means of eqs. (25), 
(18), (19) and (21). At high temperatures and sufficiently strong magnetic fields the 
ratio (m, ¥.;/ma)(t/env)? containing the coulomb collision frequency »,; bccomes 
much less than unity. Then, the major part of the heating will be produced by 
collisions between charged and neutral particles, and not by coulomb collisions, i.e., 
the frictional work of eq. (25) is reduced to 


wr aor? yao OP? 7. (50) 


This is a situation similar to that in the solar chromosphere and in cool interstellar 
clouds where the dissipation is enhanced by collisions between ions and neutral 
particles (Piddington 1954, 1956; Cowling 1956 and Lehnert 19585). 

Collisions between the reentering neutral gas and the plasma will also enhance the 
diffusion of charged particles into the loss cones of a mirror machine and reduce the 
confinement times beyond the values given by pure coulomb scattering. 

Finally, a few comments should also be made on the conduction of heat by charged 
particles. According to Kennard (1938) the heat conductivity in a gas of density n, 
mean thermal velocity 6 and a mean free path L is about 4nk@L. In the case of a 
magnetized plasma the heat conduction across a strong magnetic field will probably 
be reduced considerably, whereas it is unaltered by the field in the longitudinal 
direction (Spitzer 1956). Take an example with a hydrogen plasma with n = 1022m-3, 
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5=3 x 10’ m/s corresponding to an electron temperature of about 10’ °K and with 
Lx0.3m. Then, the heat conductivity becomes of the order of 10° watts:-m—!, 
degree~* along the magnetic field, which is greater than that of solid copper. Thus, 
at the long free paths which are actual in a discharge at very high temperatures the 
temperature gradient along the magnetic field lines is likely to be very small. These 
are about the same conditions which prevail in the exosphere, which is isothermal 
(Spitzer 1949, Bates 1954 and Nicolet 1954). 

In the next two paragraphs examples will be given of two limiting cases, one without 
heat conduction and one with infinite heat conductivity along the magnetic field. 


b. Plasma without heat conduction 


First consider a simplified situation where no thermal heat flow Q is produced by 
charged particles. Further restrict the discussion not to include elastic scattering of 
neutrals out of the confinement volume, i.e., g;,=0, 4 =0 according to sections 
II.1 and II.5.c. The ion and electron temperatures are assumed to be approximately 
equal, ie., 7, 7, = 7. According to the preceding section IT1.2.a v7 ~ v2 ~ V2  Q???. 
From eqs. (43), (41) and (46) follows that a very strong magnetic field forces the 
convection to take place only along the field lines in the rz-plane and 


d p) p) d 
qj (P/@) = (12+ ne) (p/e) * 0 7 (P/e); (51) 


where the derivative is taken along a field line. A similar expression is obtained for 
do |/dt. 

tae again a rough estimation of orders of magnitude it is found from eqs. 
(22)-(25) and from eq. (46) that the ratio E-i/en v,, H becomes very small when the 
radius of gyration is small compared to the dimensions of the confinement region, 
provided that the thermal velocity of a charged particle does not exceed the angular 
velocity Qr considerably. The bounding material walls at the ends of a field line are 
assumed to be non-conducting. Since div i=0 it is then reasonable to assume that 
the part of i in the rz-plane which is perpendicular to E also becomes small compared 
to en v,. This simply means that the ions and electrons move with only slightly 
different mass velocities in the rz-plane. 

With these approximations and from eq. (50) the energy theorem (32) becomes 


d d 
(1+f)EoQ? *—fep-BfEp~lor 7,(/e)— rig, th +N EoQ”, (52) 


where m,/m has been put equal to unity. 
Introduce the variable t=kT/[mQ?. (53) 


Eqs. (47), (45) and (52) are now reduced to 


d d e 
gr5-=—-(07) (54) 
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and (1+ fr? (5+ $f) t= (0/8) (87-9). (55) 


Two special solutions are studied. One is obtained when charge exchange is neg- 
ligible (f<1). Then eqs. (54) and (55) give the solutions 


r=r/5 (56) 
and 0/00 = (7/1); P| Po = (1/70) (57) 


If, instead charge exchange dominates over ionization (f> 1) the corresponding solu- 


tions become 
tT=2 77/3 (58) 


and e/eo= (r/o) *s — p/Po= (1/1)? - (59) 


The solutions (56) and (58) are plausible from the point of view of elementary con- 
siderations. An ion which is formed from a neutral particle by ionization or charge 
exchange will start with zero velocity and perform a cycloidal orbit in the crossed 
fields E and B. This corresponds to a mass motion of energy 4m,Q?r? and an equal 
amount of ‘“‘thermal’’ energy with two degrees of freedom. The latter is now trans- 
ferred into three degrees of freedom. Further it has also to be shared with the electron 
which is formed when an ionizing collision is concerned, but not in the case of charge 
exchange. This would give t =7?/6 and t =7?/3 for ionizing and charge exchange 
collisions, respectively. However, these values of t are too small, since a transforma- 
tion of “thermal” energy from two to three degrees of freedom requires a deflection 
of the ion by which it will “fall” to a somewhat “lower” electric potential surface 
than that “grazed” by the cycloidal orbit. This somewhat increases the values of t 
given by this elementary consideration, as is also shown by eqs. (56) and (58). 

The results obtained so far merely serve as an illustration of the energy balance and 
will not be treated further. The discussion is now continued with the case of infinite 
heat conductivity in the direction of the magnetic field, a situation which according 
to section IT.2.a is more likely to occur in plasma experiments at high temperatures. 


c. Plasma with infinite heat conductivity along the magnetic field 


The plasma atmosphere to be studied is now supposed to be isothermal along the 
magnetic field lines. This is likely to be a good approximation in a number of interest- 
ing applications. For instance, a plasma with n = 1022m-3, 7210’ °K, D=0.3 m 
and 10 percent temperature difference over a distance of 1 m along the field lines 
would have a heat flow of about 10! watts/m?. Consequently, equilibrium requires 
an almost constant temperature along the magnetic field. 

With the approximations of the preceding paragraph, except that of vanishing 
heat conductivity, eq. (47) can be integrated to 


0 = Oy exp [— (r6 — 7?) /2 ty], (60) 


where T, is given by the definition (53) and the constant temperature 7’, along the 
field line in question. Eq. (60) is analogous to the “barometer formula” where the 
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centrifugal field has been substituted for the gravitation field. An equivalent expres- 
sion has earlier been obtained by Block (1958) in connection with a study of the co- 
rotation with the sun of the interplanetary gas. 

The degree of confinement depends on the magnitude of 1, and this can be deter- 
mined from the energy balance. Since heat conduction and heat convection cannot 
be separated in this case the energy balance has to be considered for a complete 
strip like that indicated in Fig. 1, which is extended from the equatorial plane where 
r =frq up to a point r =r, where charged particles hit a material wall and get lost. 

With the notations of Fig. 1 conservation of magnetic flux requires that 


rBdl=r,B,dl, (61) 


and a volume element of the strip at r becomes dw = 22(B,/B)rjdl,ds. Then, the 
equation of continuity (12) can be integrated over the volume of the strip to give the 
total number of charged particles of either sign lost per unit time at 7; 


aN,=—r, By dy: 1y, (62) 


where index (,) indicates values at r = 7, and use has been made of relations (46) and 
(61). The bounding material wall at r=7, is assumed to be nonconductive and is 
charged negatively in such a way that ions and electrons escape at the same rate 
through the thin Debye sheet close to r =7,. This rate is given by eq. (62). The par- 
ticles escaping at the edge of the Debye sheet correspond to an energy loss per unit 
time 

dW, =2 7% By dy IT, Q? (3 t) +11). (63) 


The corresponding energy loss caused by charge exchange becomes 


dWers= 27 Bodly O | fGx +19) ds. (64) 
0 


A suitable form of expressing the energy balance is now given by eq. (34) where q 
represents the energy flow of charged particles by convection and heat conduction. 
For the strip in Fig. 1 the energy balance requires that 

| div qdw=dW, =f (w, + w, + wy) dw —adWer, (65) 
when elastic scattering of neutrals from the confinement region is again neglected. 


With the approximations already discussed in the preceding paragraphs and with 
eqs. (23), (24), (42), (44), (46) and (12) the relations 


d 
Wa © 0 (¥*V) (v°/2) sor 07 0,5 (66) 


d 
and we.¥ —pdivvy —pét+ sas (67) 
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are obtained. Combination of these relations and eqs. (50), (46), (47) and (64) gives, 
after some deductions and a partial integration where it is observed that I(0) =9, 


the following result: 
cont 68 
(6+Eh)m—O+h)AL= | Tbr —P) fds se 
0 


It should be stressed that in the deductions of eqs. (60), (68) and (56)-(59) no 
assumption has been made about the special geometry of the magnetic field or about 
the special form of the functions € and f. When 7’ =const. along a field line this 
will also be the situation for f. In the extreme case of a homogeneous magnetic field 
r =r, and for f being a constant the right hand sides of eqs. (55) and (68) vanish. 
Both equations then give the relation ; 


T=%H=N(L+f)/(5+8/). (69) 


The modulus of the right hand member of eq. (68) is less than f J, (37) + 
+|rs—77|). If the discussion is restricted to cases where ionization dominates 
over charge exchange (f<1) the result becomes 


T=n/5 (70) 
and the density and pressure distributions are 


e/eo= exp [—5 (ro—1°)/2 ri] = p/ Dy. (71) 


In cases where 7, is definitely less than 7, as in the field from a current loop, the den- 
sity o(r =1,) at the bounding wall will be reduced very strongly. The situation de- 
scribed by eq. (70) gives a thermal energy which is low enough compared to the energy 
of the rotating mass motion to secure a good plasma confinement, i.e., a thermal 
velocity which is definitely below the “escape velocity’’ in the centrifugal field. The 
difference between the results (56) and (70) is caused by the high thermal conductivity 
along the field which reduces the temperature to a value associated with the cold 
material wall at r =7,. Thereby it should be observed that, according to elementary 
kinetic theory, the heat conductivity is independent of density and is not affected 
by a very steep density gradient as that given by eq. (71). However, it should be kept 
in mind that the applicability of a fluid model becomes doubtful when the. mean 
free path is much larger than the dimensions along the magnetic field. 

An approximate expression for the energy loss by charged particles at the edge 
r =r, is obtained when the ions are assumed to enter the Debye sheet with a mean 
veolocity %;=(8k T/m)} in the longitudinal direction. This produces a particle 


flux n, %;/4 and an energy loss 


dW, =2ar,dl,-n, (4/4): (37; +i mr). (72) 


The normalized power loss from the strip counted per unit area of the equatiorial 
plane at r = ry is given by eqs. (70), (61), (72) and (71); 


: 2 
T= \/2-06b2/m)t- (By/B,) exp |-3(3-4)]- (73) 
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Also observe that, in addition to the reduction of density in the regions of the “‘holes’”’ 
of a magnetic bottle, the current loop also has the advantage of a very small geo- 
metrical cross section at the edge r = r, as shown by the factor B,/ B,. 

The results may finally be illustrated by some numerical examples. Assume a 
temperature of 7’ =2 x 10° °K. For deuterium f and h become about 0.12 and 0.07, 
respectively, and eqs. (70), (71) and (73) can be used with good approximation. With 
a density n = 3 x 10”? m- corresponding to 1 mm pressure at room temperature four 
examples have been calculated as shown in Table 1. The first corresponds to a homo- 


Table 1. The specific power loss [[, due to particle losses from a rotating plasma 
confined in a magnetic field at 7=2x10® °K and a density n=3x10" m7’. 


a (m\3 
ro/71 B,/B, 2 (=) II,/e, (watta/m®) 

Homopolar machine 

(almost homogeneous field) 1 1 1 1.2 x 104 

(mirror ratio 2) .... 1.414 0.500 0.041 4.9 x 109 
Peasteniloop. 2. let thes. i 3.03 0.054 7.45 x 10—*} 8.9 
Current loop with auxiliary field 5.50 0.115 2:30 x 10-8 2.8 x 10-22 

Se a ea | ok ee tL ae ee oes bas ae Sees Dee 


polar machine with a mirror ratio near unity, the second to a machine with a mirror 
ratio of 2, the third to a current loop where 7, is equal to twice the radius R of the 
loop and the fourth to a loop with an auxiliary field which “‘compresses” the mag- 
netic flux into a sphere (Lehnert 1959) of radius 2.5 R. In this example 7, has been 
put equal to 2.2 R. The last two examples of the table show that, as a consequ- 
ence of the large ratio r,/r, and the small cross section at 7,, the energy losses by 
escaping charged particles can be neglected in the magnetic bottle formed by a cur- 
rent loop. Other losses such as those produced by radiation will probably determine 
the energy balance. 

According to the first two examples of the same table the energy loss in an ordinary 
homopolar machine is expected to be much higher. This seems to be a plauisble 
conclusion even if the values given in the table correspond to the isotropic pressure 
distribution and should be considered as an upper limit. They are also uncertain due 
to the fact that they have been derived from a theory for small particle losses, and 
this is not the situation when 7,/7, obtains values close to unity. 

The energy input in the homopolar device described by Anderson et el. (1958) is 
of the order of 10'° watts/m? at a density of about 2 x 107? m-* and an ion energy 
of about 50 eV in deuterium (Colgate 1958). With 1/7, ~1 eq. (73) gives 7 x 10° 
watts/m? which is of the right order, even if the applicability of eq. (73) in this case 
becomes somewhat doubtful. The major part of the energy input is then likely to be 
lost by particles moving along the magnetic field to the vessel walls. The same 
conclusions can also be drawn for the modified homopolar with nonuniform field 
(Colgate 1958) where the ratio r,/7, only exceeds unity by a small amount and where 
the cross section of the vessel wall at 7 =r, is comparatively large. 


IV. The stability problem 


So far no attention has been paid to the possible existence of instabilities which 
might increase the energy dissipation catastrophically and frustrate the attempts to 
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obtain a high temperature. As far as can be judged at the present stage a sufficiently 
low plasma pressure and a sufficiently strong magnetic field is likely to give a stable 
configuration, at least as far as magnetohydrodynamic instabilities are concerned. 

In fact, the mirror machines have been found by experiments to be stable (Post 
1958). This may seem to be in contradiction with theoretical work made by a number 
of authors which have predicted the existence of ‘interchange instabilities”. How- 
ever, care is necessary when these results are applied to experiments. In a great 
number of theoretical investigations a situation has been studied where a surface 
current divides a vacuum region with magnetic field from a fieldfree plasma region. 
This is a very special problem where deformations of the bounding surface correspond 
to a displacement of all currents which produce the field discontinuity at the surface. 
A completely different situation is obtained in a plasma with a continuous pressure 
and current distribution. Here a displacement of a fluid element does not necessarily 
change the whole current distribution and the magnetic field as in the case of the 
sheet current. For the latter it has been found that an “interchange instability” 
occurs when the magnetic field lines bend concavely toward the plasma. This rule 
has been shown by Berkowitz et al. (1958) not to be applicable to mirror configurations 
with volume currents. These authors also point out that for a system with volume 
currents “instabilities can develop only because of cooperative behaviour involving 
an appreciable part of the system’. The possibility of stabilizing the system with 
external magnetic fields ‘‘anchored”’ in solid conductors also requires further inves- 
tigation. 

Consequently, it is not at all sure that a system with magnetic field lines bending 
concavely toward the plasma should be unstable. If instabilities would in fact arise 
they could possibly be suppressed by using a twisted magnetic field with “‘shear”’ as 
suggested by Spitzer (1956). The possibility of adding a toroidal field to that genera- 
ted by the current loop has earlier been discussed by Lehnert (1959). For the sta- 
tionary motion of a rotating plasma described in this paper nearly all the current 
will flow in the y-direction around the axis, provided that the losses are very small. 
The presence of an additional toroidal field b, hardly influences the forces and the 
pressure balance in such a case. If a disturbance of the axisymmetric motion arises, 
however, the shear produced by b, comes into operation. 

Finally, an enhancement of the transverse diffusion rate may also be caused by 
asymmetric and locally varying electric fields in the plasma. Experimental results 
which may indicate the existence of such a phenomenon have been observed in the 
positive column of a discharge (Lehnert 19580, d, Hoh and Lehnert 1959). On the 
other hand, no such phenomenon seems to be present in the experiments with mirror 
machines (Post 1958).'There are a number of facts which might give possible explana- 
tion to the difference between these results, and which have to be analysed further. 
The former experiments are made with a relatively cold discharge of low ionization 
degree in a narrow tube with an electric field directed along the magnetic field lines. 
The latter are made with a hot, almost fully ionized gas in a mirror machine of rela- 
tively large dimensions, where almost no longitudinal electric field is present. 


Royal Institute of Technology, Stockholm 70, December 19, 1959. 
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